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Abstract 

The  main  objective  of  this  study  is  to  develop  the  synthesis  and  optimization  of  reactor- separation  network  (RSN)  models  that  can  be  simultane¬ 
ously  solved  within  Non-Linear  Programming  (NLP)  for  a  PEM  fuel  cell  system.  The  objective  function  for  optimization  was  defined  to  minimize 
the  overall  cost  and  CO  production.  Live  alternatives  were  synthesized  to  determine  the  best  flow  chart  for  the  system  based  on  cost  and  the  output 
concentration  of  carbon  monoxide.  A  Polymer  Electrolyte  Membrane  Fuel  Cell  (PEMFC)  system  was  taken  as  the  case  study.  The  results  indicated 
that  the  optimum  specific  cost  of  a  PEMFC  stack  was  found  to  be  in  the  region  of  US$  500-700  kW-1,  while  the  specific  manufacturing  cost  and 
the  specific  investment  cost  were  calculated  at  the  range  of  US$  1000-1500  and  2500-3000  kW-1,  respectively.  Furthermore,  the  infrastructure 
investment  cost  was  determined  to  be  in  the  range  of  US$  10-30  billion  with  the  specific  cost  for  one  unit  in  the  range  of  US$  2000-4000  kW-1. 
The  results  obtained  are  comparable  with  other  studies. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Very  few  publications  exit  on  the  optimization  of  reactor- 
separation  networks  (RSN)  although  in  most  chemical  pro¬ 
cesses,  reactors  are  sequenced  by  separation  systems.  Most  of 
the  research  work  on  optimization  only  concentrates  on  the 
heat  management  [1-10]  and  reaction  system  but  the  interaction 
between  the  structure,  design  and  operation  of  both  reaction  and 
separation  units  are  not  fully  exploited  particularly  for  a  fuel  cell 
system.  The  limitation  of  this  approach  is  due  to  the  complexity 
of  the  formulation  involved  in  doing  so,  but  these  parameters  are 
very  important  in  process  system  engineering  in  order  to  ensure 
the  purity  of  the  final  product  and  to  reduce  the  overall  cost  as 
well  as  waste  generation. 

2.  Problem  statement 

The  problems  to  be  addressed  in  this  paper  can  be  stated  as 
Table  1. 
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The  main  objective  of  this  study  is  to  introduce  synthesis 
of  the  superstructure-based  approach  to  the  synthesis  of  a  RSN 
network  For  PEM  fuel  cell  system  as  follows: 

•  To  develop  reactors  with  a  continuous  reaction  and  separation 
system,  which  are  based  on  Linear  Programming  (LP)  and 
Non-Linear  Programming  (NLP)  synthesis  procedures. 

•  To  develop  a  simplified  and  user  friendly,  yet  a  general  super¬ 
structure  of  RSN  networks  in  matrix  form  and  solved  using 
MATHLAB,  and  the  overall  process  synthesis/optimization 
is  carried  out  at  different  level  of  superstructure  details  and 
model  aggregation. 

3.  Methodology 

The  optimization  problem  is  considerably  more  difficult  if 
the  variables  included  take  into  consideration  the  integer  and 
binary  (0-1)  values.  These  problems  are  referred  to  as  Mixed 
Integer  Non-linear  Programming  (MINLP)  problems  and 
usually  become  the  hardest  part  to  be  solved  unless  a  special 
structure  can  be  exploited.  However,  in  practice,  most  design 
problems  are  not  solved  using  superstructures  and  mixed-integer 
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Nomenclature 

A 

cell  active  area 

A 

linear  inequality  constraints 

Aeq 

linear  equality  constraints 

Ao 

pre-exponential  factor 

Am 

area  of  membrane 

Aa 

area  of  adsorber 

Af 

active  area  of  fuel  cell 

ATR,  Rl 

auto-thermal  reformer 

c/Cf 

degree  of  purification 

f*eq 

non-linear  equality  constraints 

c 

non-linear  inequality  constraints 

CA 

methanol  cost 

Cb 

water  supply  cost 

Cc 

oxygen  supply  cost 

Cr-  1 

cost  for  Reactor  1 

Cr-2 

cost  for  Reactor  2 

Cr- 3 

cost  for  Reactor  3 

Cr- 4 

cost  for  Reactor  4 

Cs  - 1 

cost  for  unit  separation  1 

Cs  —  2 

cost  for  unit  separation  2 

Cm ,  Cs,  Co2,  Cco  molar  concentration  for  methanol, 
steam,  O2,  CO 

^Utilities 

cost  for  utilities 

CT 

total  investment  cost 

D 

diameter 

Ea 

activation  energy 

F 

Faraday  constant  (96,487  C) 

fix) 

objective  function 

Fmo 

feed  concentration  of  methanol 

fp 

cost  factor  for  pressure 

fm 

cost  factor  for  material 

FC,  R5 

fuel  cell  stack 

g(x) 

equality  constraints 

h(x ) 

inequality  constraints 

[H] 

Hessian  matrix 

i 

current 

io 

current  density 

^Cost 

cost  index 

J 

activity  of  inequality  constraints 

k\,  k2,  kM,  ks,  ko2Jt\\,  kco  reaction  constant  reaction  1, 

2,  methanol,  steam,  oxygen,  hydrogen,  carbon 
monoxide 

L 

length 

Lb,  Xlow 

lower  bound  limit 

m 

number  of  equations 

Mo,  Ms 

ratio  of  methanol  to  hydrogen  and  steam 

MR,  R4 

membrane  reactor 

n 

number  of  variables 

^fuel  cell 

fuel  cell  efficiency 

P 

pressure 

Pm 

partial  pressure  for  methanol 

pw. 

partial  pressure  for  hydrogen 

POX 

partial  oxidation 

PROX,  R3  preferential  oxidation 
PSA,  S2  pressure  swing  adsorber 
Qu  Q2,  63  heat  balance  for  reactors  1-3 
ri,  r2,  rM,  rs ,  ro2,  Hd  rco  rate  of  reaction  for  reaction 
1 ,  2,  methanol,  steam,  oxygen,  hydrogen,  carbon 
monoxide 

Sh2  hydrogen  selectivity 

S  search  direction 

S\,  S2,  S3,  S4  selectivity  for  reactors  1-4 

SR  steam  reforming 

To  initial  temperature 

Tu  T2,  T3  temperature 

TCM,  S 1  tubular  ceramic  membrane 

Up,  Xup  upper  bound  limit 

Vceii  voltage  of  single  cell  (V) 

V,  Vi,V2,V3,  Vatr,  Vwgs,  Uprox  volume  of  reactor 
w  width 

WGS,  R2  water  gas  shift  reactor 

Xfj  mole  fraction  of  component  i  in  feed  stream 

x\,  x2,  . . xn  design  parameters 

Xa  conversion  of  component  A 

Greek  letters 
a  step  calculation 

a  ideal  separation  factor 

£w  volume  expansion  factor 

0  pore  diameter 

X  the  conversion  of  H2  to  H2O  in  PROX 

0  degree  of  separation 

p  density 


optimization  algorithms  only  [1 1] .  In  addition,  the  mixed  integer 
formulation  of  the  optimization  problem,  in  principle,  permits 
the  optimizer  to  select  simultaneously  the  best  flow  sheet  and 
optimize  it  with  respect  to  its  continuous  variables  [12,13].  Due 
to  that,  in  this  study  the  initial  stage  generations  of  the  RSN 


Table  1 

Problem  statements 


Given 

Find 

Rate  of  chemical  reaction  in  the 

The  reactor  network  (i.e.  size  of 

reactors 

reactors,  feeding,  conversion  of 

Feed  flow  rates  and  composition  of 

reactants,  etc.) 

The  separation  network  strategy 

feed  streams 

(e.g.  the  size  of  the  membrane 

Degree  of  separation  for  separators, 

and  adsorber,  and  the  sequence) 
The  interconnection  between  the 

i.e.  membrane  and  adsorber 

two  networks 

Physical  and  chemical  properties  of 

Selecting  the  best  alternative  by 

components 

synthesizing  the  RSN  based  on 

A  performance  criterion  (e.g.  total 
cost,  yield,  selectivity  of  desired 
product) 

Mass  balance  of  the  system 

Energy  balance  of  the  system 

minimum  cost  and  product  purity 
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networks  are  simplified  to  the  Non-Linear  Programming  (NLP) 
procedures. 

Although  the  NLP  approach  is  computationally  more  exten¬ 
sive,  it  has  been  proven  to  be  a  more  robust  and  reli¬ 
able  method  [14-16].  The  NLP  problems  in  the  study  were 
solved  using  optimization  tools  in  MATHLAB.  MATHLAB 
is  a  standard  tool  introduced  for  technical  computing  and 
was  originally  developed  as  a  matrix  laboratory.  The  major 
element  of  MATHLAB  is  a  matrix.  Mathematical  manipu¬ 
lations  of  matrices  are  very  conveniently  implemented  and 
easy  to-use  as  an  interactive  environment.  To  these  matrix 
manipulations,  MATHLAB  has  added  a  variety  of  numeri¬ 
cal  method  functions,  symbolic  computations  and  visualization 
tools  [30,31]. 

SQP  has  emerged  as  a  very  popular  and  powerful  method 
for  solving  NLP  problems  in  process  optimization  [12,15,17]. 
Fig.  1  summarizes  the  flow  chart  for  the  development  of  an  NLP 
algorithm  for  SQP.  [H]  or  Hessian  matrix  is  referred  to  as  A2f(Xi) 
and  S  is  referred  to  as  search  direction. 

Optimization  techniques  are  used  to  find  a  set  of  design 
parameters,  i.e.  x={x\,  X2,  . . .,  xn},  that  can,  in  some  ways, 
be  defined  as  optimal.  The  development  of  algorithms  for  syn¬ 
thesis  and  optimization  of  the  RSN  in  this  study  were  based  on 
the  modified  finite  difference  derivatives  in  sparse  matrix.  These 
algorithms  can  be  used  to  solve  function  with  multiple  variables 


Fig.  1.  Flowchart  of  NLP  Algorithm  for  SQP. 


with  constraints  as  follows: 


Minimize  f{x)  (1) 

s.t. 

c(x)  =  0  (2) 

Ceq(x)  =  0  (3) 

Ax  <  b  (4) 

AcqX  <  bQ q  (5) 

lb  <  x  <  ub  (6) 


with  v,  b ,  Z?e q,  lb  and  ub  are  vector  in  the  form  of  sparse  matrix, 
A  and  Aeq  are  the  sparse  matrix  returns  the  values  for  linear  con¬ 
straints,  c(x)  and  ceq(v)  are  the  function  for  non-linear  constraints 
translates  the  vector  and  the/(v)  is  the  function  that  translates  in 
scalar. /(x),  c(x)  and  ceq(v)  are  non-linear  function. 

4.  System  description 

Fig.  2  shows  the  schematic  diagram  of  a  PEMFC  system. 
The  main  electric  power  generator  in  the  system  is  the  fuel  cell 
stack  (that  contains  the  anode,  cathode,  PEM  and  cooling  plates) 
as  shown  in  the  diagram.  Fuel  gas  containing  H2  is  fed  to  the 
anode  side  and  oxygen  (or  air)  is  fed  to  the  cathode  side  of 
the  stack.  The  spent  fuel  gas  and  the  gas  leaving  the  fuel  cell 
stack  are  referred  to  as  the  anode  exhaust  and  cathode  exhaust, 
respectively.  For  convenience  in  describing  the  fuel  cell  system 
configuration  as  shown  in  Fig.  2,  the  following  paragraph  and 
Table  2  shah  explain  the  layout  of  the  system  in  terms  of  the 
flows  of  oxygen  (or  air),  fuel,  processed  water,  exhaust  gases 
and  stack  coolant. 

Beginning  at  the  right-bottom,  ambient  air  is  compressed  to 
3  atm  and  passes  through  an  oxygen  concentrator  to  separate 
the  O2  from  N2.  This  stream  will  split  into  two  main  streams. 
The  smaller  oxygen  (or  air)  flow  is  preheated  and  fed  to  an  ATR 
reactor.  On  the  other  hand,  the  larger  oxygen  (or  air),  which  flows 
to  the  HUM  1 ,  is  humidified  and  cooled  by  injecting  liquid  water 
prior  to  being  supplied  to  the  cathode  inlet  of  the  fuel  cell  stack. 
The  fuel  flow  circuits  begin  at  the  fuel  tank  (methanol).  From 
the  tank,  it  is  pumped  to  the  fuel  vaporizer  where  it  is  heated 
and  vaporized  by  the  hot  stream  exhaust  gas.  It  is  then  fed  to 
the  ATR  where  it  reacts  with  the  preheated  air  and  superheated 
steam  to  yield  the  raw  reformate.  The  hot  reformate  is  used  to 
superheat  the  ATR  steam  feed.  It  is  then  mixed  with  additional 
water  and  fed  to  a  WGS  reactor.  The  fuel  gas  leaving  this  reactor 
will  pass  through  the  water  adsorber  unit.  The  purpose  of  this 
unit  is  to  remove  the  excess  water  in  the  stream  before  entering 
the  membrane  units.  The  flow  that  leaves  the  membrane  unit  is 
then  cooled  by  lowering  the  temperature  of  the  inlet  oxygen  (or 
air)  and  subsequently  vaporized  the  feed  methanol  to  the  ATR 
before  entering  the  adsorber.  Finally,  the  purified  hydrogen  fuel 
flow  will  be  humidified  to  100%  by  a  HUM  2  and  will  enter  the 
anode  side  of  the  stack. 

The  exhaust  gas  circuit  is  best  traced  from  the  upper  left 
corner  in  Fig.  2.  The  fuel  cell  cathode  exhaust  is  cooled  to  recover 
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Coolerl 


Fig.  2.  The  proposed  schematic  diagram  of  PEMFC  system  in  this  study. 


Table  2 

System  layout  and  alternatives 


Items 


Description 


Reaction 


Auto -thermal  reformer  (ATR) 


Producing  H2  from  methanol  (A)  and  steam  (B) 
with  co-feeding  with  oxygen  (C) 


Water  gas  shift  reactor  (WGS) 

Preferential  oxidation  (PROX) 
Tubular  ceramic  membrane  (TCM) 
and  pressure  swing  adsorber  (PSA) 
Fuel  cell  stack  (FC) 

Membrane  reactor  (MR) 


Secondary  hydrogen  producer  and  primary  CO 
clean-up  system 

Alternative  to  WGS  in  removing  CO 
Operated  in  parallel  to  yield  the  product  purity  as 
99.9%  hydrogen  and  CO  less  than  20  ppm 
Consists  of  three  major  components;  anode,  cathode 
and  platinum-catalyzed 
Combines  ATR,  WGS  and  TCM 


Endothermic  steam  reforming:  CH3OH  +  H2O  CO2  +  3H2 
exothermic  partial  oxidation, 

CH3OH  +  0.5O2  -*  C02+2H2+(CO) 

C0  +  H20^C02+H2 

CO  +  o2  ->  co2 


Anode:  2H2  — >  4e  +  4H+ ;  cathode:  4e  +  4H+  +02^  2H2O 
Overall:  2H2  +  02^2H20 


water  and  purged  to  the  air.  The  last  subsystem  is  the  stack 
coolant  circuit  that  is  used  to  remove  the  waste  heat  from  the 
fuel  stack.  Coolant  from  the  tank  is  pumped  to  the  cooling  plates 
in  the  fuel  cell  stack.  The  heated  coolant  then  flows  to  the  main 
radiator  and  then  back  to  the  coolant  tank. 


5.  Model  development  of  RSN  network 

Fig.  3  shows  the  alternative  configurations  of  the  process  units 
in  the  system  in  block  diagram  superstructure  while  Fig.  4  shows 
the  superstructure  embedded  with  the  five  alternatives. 


(a)  Alternative  1  Alternative  2 


Alternative  3  Alternative  5 

Alternative  4 


Fig.  3.  Four  structural  alternatives  configuration  of  the  units  in  the  system. 
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The  following  paragraph  presents  the  mathematical  formula¬ 
tion  for  RSN  network  by  defining  the  variables  for  the  system: 

i.  The  existence  of  each  unit  of  reactor  and  recycle  streams 
(e.g.  for  NLP  problem,  the  y  variable  is  set  to  1  or  0  for  the 
continuous  reaction). 

ii.  The  existence  of  the  separation  units  (e.g.  for  NLP  problem, 
the  z  variable  is  set  to  1  or  0  for  the  separation  units). 

iii.  The  stream  flow  rate  (i.e.  F\ ,  F 2,  7%  . . .,  F35). 

iv.  The  component  compositions  (molar  fractions)  at  the  inlet 
of  each  reactor  unit  (i.e.  X4 x^j,  xsj,  i=A-C )  and  at  the 
outlet  of  each  reactor  unit  (i.e.  x^j,  x-/j,  xgj,  i—A-C ). 

v.  The  volume  of  the  reactor  units  (i.e.  Vrj,  Vr?2,  Vr,3). 

vi.  The  feed  flow  rate  of  components  A-C  and  the  final  flow 
rate  of  components,  D-F. 

vii.  The  component  composition  at  the  inlet  of  TCM  (i.e.  X22 ,u 
i  =  B,D,E,F). 

The  set  of  constraints  consists  of 

i.  mass  balances  for  splitters  of  feed  and  reactor  units; 

ii.  component  balances  for  mixers  at  the  inlets  of  the  reactor 
units; 

iii.  component  balances  for  mixers  at  the  inlet  of  TCM  and 
PSA; 

iv.  component  balances  at  the  mixers  prior  to  recycles  alter¬ 
natives; 

v.  component  balances  around  each  reactor; 

vi.  component  balance  for  each  separation  units; 

vii.  summation  of  mole  fractions; 

viii.  logical  constraints  and  volume  constraints; 

ix.  non-negativity  and  integral  constraints. 


The  reactions  involved  in  the  reactor  systems  are  as  the  fol¬ 


lowing  mechanism: 

Reactor  1  :  A  +  B  +  C  — ->  D  +  E  +  F  (7) 

Reactor  2  :  F  +  B-^D  +  E  (8) 

Reactor  3  :  E  +  C  — F  (9) 

and 

d  +  c^Ib  (10) 

Reactor  4  :  A  +  B  +  C— — *D  +  E  +  F  (11) 

and 

F  +  B^Id  +  E  (12) 

Reactor  5  :  D  +  C-L  B  (13) 


where  D  is  the  desired  product,  A-C  are  the  feed  and  E  and  F 
are  the  by-product. 

Table  3  shows  the  kinetic  constants  for  the  reactions  involved, 
while  Table  4  presents  the  summary  of  the  shortcut  design 
method  models  used  for  the  design  of  various  units  in  the  system. 
More  detailed  information  can  be  found  in  the  cited  references. 
Besides  that,  Table  5  summarizes  the  operating  conditions  for 
each  unit  in  the  system,  meanwhile  Table  6  shows  the  cost  param¬ 
eter  for  the  components  in  the  system. 

5.7.  The  objective  function 

From  the  overall  reactions,  D  (H2)  is  the  desired  product,  A 
(methanol),  B  (water/steam),  C  (O2)  are  the  feed  and  E  (CO2) 
and  F  (CO)  are  the  by-products.  The  objective  is  to  minimize 
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Table  3 

Activation  energy,  pre-exponential  factors  and  reaction  constants 


Table  6 

Cost  parameter  [27-29] 


Equations 


Values 


Components 


Cost  (US$)/unit 


r  _  u  rl.l5  O-Xm)(Ms-Xm)°-5(M0-Xm)0-25 
rATR  -  hCM0  - w+£.„Zm)..75 - 


— AVGS  =klC 


2 

COO 


(l-Xco)(M's-Xco) 
(i+sw  Xco)2 


kC0k02PC0PQ^ 

r3  0.5.kCopco+ko2po2 


r  4 


(Kco/1+KdKCOPm/P^2  )( Kdpm/p -pco) 

a 


ATR  [19] 

•Ea  =  100.9  kJ  mol-1 
•A0  =  1.9  x  1012 
mol(min  gcat  kPa)  - 1 

WGS  [20] 

•Ea  =  50  kJ  mol-1 
•A0  =  1.3  x  107 
mol(min  gcat  kPa)  - 1 

PROX  [21] 

•kco 

•  Ea  =  52.2  kJ  mol-1 

•  A0  =  8.6  x  108 
mol(min  gcat  kPa)  - 1 

•ko2 

•Ea  =  6.02  kJ  mol- 1 
•A0  =  1.75  x  107 
mol(min  gcat  kPa)  - 1 

Membrane  reactor  [22] 

•^aCO  =  V.7x  10-10 

mol  cm-3  satm 
•Kqo  =5.77  x  10-3  atm-1 
•Kd  =  76.3  atm2 


Table  4 


Shortcut  design  models 


Units 

Models 

ATR  [23] 

Plug  flow  reactor:  adiabatic  and  non-isothermal 

WGS  [23] 

Plug  flow  reactor:  isothermal 

PROX  [21] 

Plug  flow  reactor:  isothermal 

TCM  [24] 

The  length  of  the  membrane  is  simplified  to  number 
of  unit  (NTU)  and  height  of  transfer  unit  (HTU) 

PSA  [25] 

The  PSA  is  simplified  by  using  the  Klinkenberg 
solution  and  introducing  Daud  bed  utilization  factor 

Stack  [26] 

Mass  balance  and  the  polarization  parameter 

Table  5 


System  operating  condition  for  base-case  scenario 


Parameter 

Values 

Reformer  temperature 

250  °C 

Reformer  pressure 

2  bar 

Feed  flowrate  of  MeOH 

4  mol  min- 1 

MeOH:  steam 

1:1.3 

MeOH:02 

1:0.25 

WGS  temperature 

250  °C 

WGS  pressure 

5  bar 

CO:  steam 

1:30 

TC  membrane  temperature 

250  °C 

TC  membrane  H2  outlet  temperature 

1  bar 

PSA  temperature 

70  °C 

PSA  pressure 

6  bar 

Stack  pressure 

2  bar 

Stack  temperature 

70  °C 

H2/02  ratio  in  the  stack 

1.5 

Methanol 

ATR 

WGS 

PROX 

TCM 

Pd  MR 

PSA 

Water  adsorber 

Air  purifier 

Compressor 

Heat  exchanger 

Pump 

Stack 

0.08  g  mol  s  1 

29.5  x  V^t6r  x  fp  x  fm  x  7kos 

29.5  x  x  fp  x  fm  x  7kos 

29.5  X  UpRQx  x  fp  x  fm  x  Aos 
(0.93/ n0-7)  x  m  x  fp  x  7kos 

1000  x  V0  6  x  7kos 

1.54  x  10-3  x  Lom  x  Dl05  x  fp  x  7kos 

1.54  x  10-3  x  L0-81  x  Dl05  x  fp  x  7kos 

1.68  x  10-3  x  L0-78  x  D0-98  x  fp  x  7kos 

18.5  X  P°-77  X  /jenis  X  7kos 

3.3  x  A0-0024  x  7kos 

109 

0.35  m-2 

Investment  cost  (US $  x  103) 

ATR 

93.4  x  Vi 

WGS 

145  x  V2 

PROX 

160  x  V3 

TCM 

2116.8  x  feed  flow  rate 

PSA 

39  x  feed  flow  rate 

Air  purifier 

39  x  feed  flow  rate 

Water  adsorber 

39  x  feed  flow  rate 

Compressor 

58.6  x  Pm 

Stack 

440  x  W 

the  cost  and  the  production  of  E.  The  total  investment  cost,  Ct 
is  given  as  below: 

Total  cost, 

Cj  =  Ca  +  Cb  +  Cc  +  Cr- i  +  Cr—2  +  Cr- 3  +  Cr-4 


+  Cs-1  +  Cs-  2  +  ^Utilities  (14) 

CA  =  Fix  0.08  x  9.56  x  106  (15) 

CB  =  8.0  x  106  x  6.3  x  (F2  +  Fl0),  (16) 

Cc  =  3  x  13  x  103(F3  +  Fi 7  +  F3 0),  (17) 

CR-i  =  93.4  x  103  x  Vi  +  Qi,  (18) 

CR- 2  =  145  x  103  x  V2  +  Qi,  (19) 

Cr- 3  =  160  X  103  X  y3  +  <23,  (20) 

Cr- 4  =  F29X29,D  x  9.5  X  106  (21) 

Cs- 1  =  F22  X  21 16.8  X  103  (22) 

Cs-2  =  F26  X  39  X  103  (23) 


Cutilities  =  [39  x  103(/7i2  +  Fi 9  +  F32)\ 

+  [0.175  x  106  x  8  x  330]  (24) 


Vi  = 


V2  = 


F\ 


1  -  e)fci  jq 


‘ x 1  (1  +evXi)dX' 
CAo(l-Xi)  ; 

(1  +ewX2)dX 


\l/k2p(l  -  s)CFo,s  Jo 


(25) 

(26) 


(1  -  X2) 
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v  =  +16*16+  fXl  (l  +  ewX3)2dX 

3  x/skipil  —  s)Cp0  l6  Jo  (1  -  X3)(A  -  X3)’ 

01  =  F5x5,d0.2)  +  F5x 5, £(-291)  -  [Fi(-201) 

+  F2(1.3)  +  F3(-233)],  (28) 

02  =  F\\xu,d{— 359)  +  Fiixh,£(— 384) 

-  [F8x8,H-101)  +  Fio(-231)],  (29) 

03  =  Fi8xi8,B(-12.9)  +  Fi8*i8+(-16.8) 

-  [Fiv(— 1.1)  +  F16xi6,d(-10.3)  +  Fi6x16,£(-10.6)] 

(30) 


Alternatives 


Fig.  5.  Comparison  of  specific  manufacturing  and  investment. 


[Fi  x  7.6  x  105]  +  [5  x  107(F2  +  F10)]  +  [3.9  x  104(F3  +  FX1  +  F30)]  +  [9.4  x  104Vi  +  96.940i] 
+  [1.2  x  105V2  + 96.9402]  + [1.2  x  105  V3  +  25O03]  +  [7  x  108F29jt29,D]  +  [2  x  106F22] 

+  [4  x  104F26]  +  [4.0  x  104(F12  +  FW  +  F32)]  +  [6  x  109] 


(31) 


The  mathematical  model  shown  above  is  an  MINLP  problem 
since  it  has  both  continuous  and  binary  variables  and  non-linear 
objectives  and  constraints.  The  binary  variables  participate  lin¬ 
early  in  the  objective  and  logical  constraints.  However,  due  to 


[g(X)]m  <  o 

xlow  <  x  <  xup 

With  that,  the  objective  function  is: 


Min/(X)  = 


7.6  x  105xi  +  5  x  107x2  +  5  x  107xi4  +  3.9  x  104x3 
+3.9  x  104jc24  +  3.9  x  104jc45  +  9.4  x  104Vi  +  96.9401 
+  1.2  x  105V2  +  96.94  02  +  1.2  x  105  V3  +  25O02 
+7  X  10S3V61 JC98  +  2  X  106X57  +  4  X  104X6o  +4.0  X  1043C19 
+4.0  x  104a:29  +  4.0  x  104x46  +  6  x  109 


(34) 

(35) 


(36) 


the  complexity  of  the  MINLP,  the  problem  is  solved  using  the 
NLP,  where  it  considers  only  the  continuous  variables  while  the 
binary  variables  {y,  z }  are  set  to  0  or  1.  Constraints  (i),  (vii)  and 
(viii)  are  linear  while  the  remaining  constraints  are  non-linear. 
The  non-linearity  in  constraints  and  objectives  function  is  bilin¬ 
ear  type. 

Design  parameters: 

c 

k  1 , ^2, ^3 7?fuelcell >  A  A  ^  52  ,  53  ,  54,  A,  0 ,  — , 

Cp 

V,  In ,  I,  0,  P,  Ti,  72,  73,  P\,  P2,  P3,  w,  L 
Design  variables: 

P\,  ^2,  •••  ^33  >  X\i ,  X2J  .  .  •  %33,i  ?  V\,  V2,  V3,  Q 1,  Q  2,  Q  3, 

Am,  Aa,  Ap, 
where  i-A-F. 

Associating  the  design  variable  x\-x%z  with  F\,  F2,  . . .  F33, 
*l,i,  *2,i  •  •  •  X33J,  vu  V2,  V3,  Q 1,  Q2,  03,  Am,  Aa,  Af.  Then 
the  problem  can  be  expressed  as  the  following: 

Minimize  f{X),  [ X]n  (32) 

s.t.: 

(33) 


6.  Result  and  discussion 

The  model  development  for  the  synthesis  and  optimization 
processes  in  this  study  was  based  on  the  power  output,  P. 

6.1.  Process  synthesis 

Five  chosen  alternatives  are  used  in  the  synthesis  process,  as 
explained  in  Section  5.  The  summarized  results  for  the  synthesis 
and  optimization  processes  are  as  given  in  Table  6. 

Figs.  5-7  are  plotted  based  on  the  results  from  Table  7. 
Figs.  5  and  6  show  that  the  alternative  2  gives  the  highest  specific 
manufacturing  and  investment  costs  up  to  US$  4000  and  1900, 


1  2  3  4  5 

Alternatives 

Fig.  6.  The  comparison  of  infrastructural  investment  capital. 


[h(X)],  =  o 
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Alternatives 

Fig.  7.  The  comparison  of  output  CO  concentration. 

respectively,  whereas  the  alternatives  1,  3,  4  and  5  give  lower 
costs  in  the  range  of  US$  2600-3000  and  1200-1500  kW-1, 
respectively.  This  is  because  of  the  larger  cost  of  the  PROX 
reactor  obtained  compared  to  those  of  the  ATR  and  the  WGS 
reactors  as  the  catalyst  for  the  PROX  reaction  is  very  selective 
and  expensive. 

In  Fig.  6,  it  is  observed  that  the  cost  calculations  for  the  infras¬ 
tructural  investment  for  alternative  2  comes  up  with  a  figure 
amounting  to  US$  48  billion  while  the  figures  for  alternatives  1, 
3, 4  and  5  are  only  in  the  range  of  US$  10-30  billion  (Fig.  6).  As 
the  main  objective  of  the  synthesis  is  to  minimize  its  investment 
cost  as  well  as  the  carbon  monoxide  production,  the  next  selec¬ 
tion  has  to  be  based  on  the  CO  output.  From  Fig.  7,  alternative 
1  gives  the  lowest  CO  output,  i.e.  12.6  ppm,  followed  by  alter¬ 
native  4  (included  a  membrane  reactor  unit),  i.e.  13.5  ppm,  and 
alternative  3  of  16.5  ppm.  On  the  contrary,  alternative  5  gives 
zero  CO  output  because  only  hydrogen  and  oxygen  are  used  as 
its  feed  without  any  hydrocarbon  component.  Besides  the  high¬ 
est  cost  that  the  alternative  2  incurs,  it  also  produces  the  highest 
CO  content  of  43.1  ppm,  causing  it  to  be  eliminated  from  the 
selection. 

Therefore,  alternative  1  is  the  best  choice  because  it  gives 
lower  CO  output  and  specific  manufacturing  and  investment 
costs  compared  to  the  other  three  alternatives.  Although  alter¬ 
native  5  gives  a  zero  CO  output  but  its  specific  manufacturing 
and  investment  costs  as  well  its  infrastructural  investment  cost 
are  too  high,  i.e.  US$  1350,  2950  and  28  billion,  respectively, 
compared  to  the  ones  incurred  by  alternative  1,  i.e.  US$  1200, 
2600  and  10  billion,  respectively. 


Table  8 

Input  and  output  parameter  for  5  kW  power  output  (alternative  1) 


Parameters 

Values 

Net  power  output  (kW) 

2.38 

System  efficiency  (%) 

47.57 

Methanol  input  (lmin-1) 

87.3 

Water  input  (1  min- 1 ) 

113 

Oxygen  input  (1  min  - 1 ) 

21.8 

Concentration  of  CO  output  (ppm) 

12.6 

6.2.  Optimization 

Table  8  presents  the  summarized  results  after  optimization 
for  alternative  1  with  power  output  of  5  kW,  with  reference  to 
the  operation  parameters  in  Table  5,  while  Tables  9  and  10  give 
the  detailed  design  and  economic  parameters. 

From  the  synthesis  and  optimization  processes,  the  following 
has  been  identified  in  this  study: 

•  The  WGS  reaction  is  required  in  a  fuel  cell  system  as  the 
primary  CO  remover,  followed  by  the  membrane  or  PROX 
reaction  that  acts  as  the  secondary  CO  remover  along  with  an 
adsorbent  system.  This  finding  repudiates  earlier  statements 
from  [19,30],  which  claim  that  the  WGS  reactor  is  not  required 
as  one  of  the  CO  removal  systems  in  the  fuel  processing  unit 
because  the  main  problem  with  the  WGS  reactor  is  the  low 
operating  temperature  and  the  large  reactor  size.  Instead,  it 
agrees  with  the  study  put  forth  by  refs.  [31-34],  which  prove 
that  the  WGS  reaction  is  an  integral  part  of  the  PEMFC  sys¬ 
tem  as  a  CO  removal  unit  and  a  secondary  source  of  hydrogen. 
Apart  from  that,  this  study  has  also  identified  that  the  oper¬ 
ating  pressure  can  be  used  to  control  the  size  of  the  WGS 
reactor  [23]. 

•  This  study  has  also  proven  that  the  combined  WGS- 
membrane  reaction  works  better  and  is  more  cost-effective 
compared  to  the  one  by  the  PROX-membrane  combination. 
This  is  because  the  reaction  and  effectiveness  of  the  PROX 
reactor  are  very  much  dependent  on  the  choice  of  catalyst 
with  a  high  CO  oxidation  selectivity  and  this  type  of  cat¬ 
alyst  is  expensive.  The  catalyst  for  the  PROX  reaction  has 
characteristics  that  prevent  H2  selection  by  CO  oxidation  but 
H2  oxidation  will  take  over  when  the  partial  pressure  of  CO 


Table  7 


Summarized  result  from  process  synthesis 


Parameters 

1 

2 

3 

4 

5 

Power  output  (kW) 

5 

5 

5 

5 

5 

Net  power  output  (kW) 

2.38 

2.38 

2.38 

2.38 

3.25 

System  efficiency  (%) 

47.57 

47.57 

47.57 

47.57 

65.0 

Methanol  input  (1  min- 1 ) 

87.3 

92.3 

35.5 

96.7 

0 

Water  input  (1  min- 1 ) 

113 

120 

46.1 

145 

0 

Oxygen  input  (lmin-1) 

21.8 

23.1 

8.87 

32.9 

125 

CO  output  (ppm) 

12.6 

43.1 

16.5 

13.5 

0 

Working  capital  (US$) 

1980 

3108 

1956 

2000 

1516 

Specific  manufacturing  cost  (US$  kW-1) 

1210 

1900 

1194 

1220 

1347 

Specific  investment  cost  (US$  kW-1) 

2639.6 

4144.1 

2608.2 

2665.7 

2942.2 

Infrastructural  investment  capital  (billion  US  $) 

10.3 

47.7 

10.5 

10.7 

28 
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Table  9 


Detailed  design  parameters  (alternative  1) 


Parameters 

Value 

ATR 

Pressure  (bar) 

1 

Temperature  (°C) 

250 

Conversion  (%) 

98.7 

Reactor  volume  (cm3) 

148 

Catalyst  density  (m3  kg-1) 

1.1 

Voidage 

0.5 

Rate  constant  (mol  kg-1  s-1) 

0.038 

WGS 

Pressure  (bar) 

5 

Temperature  (°C) 

220 

Conversion  (%) 

90.6 

Reactor  volume  (cm3) 

325 

Catalyst  density  (m3  kg-1) 

1.1 

Voidage 

0.5 

Rate  constant  (mol  kg-1  s-1) 

0.091 

TCM 

Inlet  pressure,  PM1  (bar) 

5 

Outlet  pressure,  PM2  (bar) 

1 

r  =  PM2/PMl 

0.2 

Temperature  (°C) 

250 

Degree  of  separation 

0.3 

Pore  diameter  (nm) 

1.5 

Membrane  thickness  (m) 

0.0005 

Membrane  module  diameter  (m) 

0.15 

Membrane  length  (m) 

0.36 

PSA 

Pressure  (bar) 

6 

Temperature  (°C) 

25 

Separation  factor  (c/cp) 

0.01 

Daud  utilization  factor  (5) 

0.2 

Superstantial  velocity  (cms-1) 

5 

Equilibrium  constant 

6500 

Mass  transfer  constant  (cm-1) 

0.02 

Bed  diameter  (m) 

0.10 

Bed  length  (m) 

0.177 

Breakthrough  time  (s) 

3.3 

FCS 

Single  cell  voltage  (V) 

0.7 

Current  (A) 

200 

Number  of  cell 

36 

Current  density  (A  cm-2) 

0.8 

Cell  active  area  (cm2  per  cell) 

250 

Length  (cm) 

15.8 

Width  (cm) 

15.8 

decreases.  In  practice,  the  reaction  of  PROX  selection  for  the 
CO  oxidation  will  be  more  difficult  with  the  decrease  in  CO 
concentration  [19,21,35].  In  addition,  the  combined  PROX- 
membrane  reaction  will  incur  a  higher  cost,  i.e.  two  times 
more,  compared  to  the  one  by  the  WGS-membrane  combina¬ 
tion. 

•  From  this  study,  the  pressurized  oscillation  adsorbent  system, 
PSA,  is  found  to  be  essential  in  CO  removal  before  the  flow 
enters  the  fuel  cell  stack.  This  is  because  the  PSA  is  very  effi¬ 
cient  in  terms  of  product  purity,  product  cost,  operating  cost 
and  capital  cost  as  proven  by  refs.  [36,37].  Even  though  the 
PSA  design  is  complicated  with  large  number  of  parameters 


Table  10 


Economic  parameters 


Components 

Cost  (US$) 

ATR 

133.9 

WGS 

214.8 

Compressor 

370.0 

Fuel  feeding  pump 

106.0 

Recycle  pump 

106.0 

TCM 

372.2 

PSA 

16.8 

Air  purifier 

6.9 

Water  separator 

10.5 

Heat  exchanger  unit 

389.0 

PEMFC  stack 

2678.6 

Manufacturing  cost 

6043.1 

Installation  cost 

9064.7 

Coat  for  1  unit  of  PEM  fuel 

13198.1 

Specific  manufacturing  cost 

1208.6  kW"1 

Specific  investment  cost 

2639.6  kW“‘ 

Infrastructural  investment  capital 

10.3  billion 

Electricity  cost 

0.04  kWh 

Table  1 1 

Comparison  Of  feed  flow  rate  of  methanol  with  other  studies 

References 

Power  output  (P,  kW) 

Value  (1  h  1 ) 

[30] 

5 

1.44 

[23,26]  (before  optimization) 

5 

1.86 

This  study  (after  optimization) 

5 

1.42 

that  must  be  considered  throughout,  this  has  been  successfully 
overcome  by  introducing  a  simple,  short  cut  design  method 
for  the  PSA  [23]. 

•  The  membrane  reactor  system  (alternative  4)  is  also  a  very 
suitable  consideration  for  the  fuel  cell  system  because  it  gives 
highly  similar  performance  and  cost  comparisons  to  alterna¬ 
tive  1 .  In  addition,  it  can  produce  smaller,  more  compact,  and 
lighter  fuel  cell  system  as  suggested  by  refs.  [22,38,39]. 

6.3.  Model  verification 

In  general,  validation  for  the  models  developed  for  process 
synthesis  and  optimization  is  carried  out  based  on  other  studies 
and  the  model  development  with  Microsoft  Excel  before  opti¬ 
mization  [23,26]. 

Table  1 1  gives  the  fuel  feed  value  required  by  the  fuel  cell 
system  with  5  kW  of  power  output.  It  is  found  that  this  value  is 

Table  12 

The  comparison  of  system  efficiency  with  other  studies 
References  Efficiency  (%) 


[40] 

32-40 

[41] 

40-60 

[42] 

40 

[43] 

40-60 

[44] 

40-60 

[26]  (before  optimization) 

45.5 

This  study  (after  optimization) 

47.5 
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Table  13 


The  comparison  of  specific  stack  cost  with  other  studies 


References 

System 

Cost  per  kW  (US$) 

[42] 

PEMFC 

500-800 

[48] 

PEMFC 

500-1000 

[43] 

PEMFC 

400-1000 

[49] 

PEMFC 

750-1000 

[45] 

PEMFC 

200-300 

[46] 

PEMFC 

400-700 

This  study 

PEMFC 

500-700 

Table  14 

The  comparison  of  specific  manufacturing  cost  with  other  studies 

References 

System 

Cost  per  kW  (US$) 

[40] 

PEMFC  (light  industry) 

1000-2000 

[48] 

PEMFC  (vehicles) 

1000-2000 

This  study 

PEMFC  (vehicles) 

1000-1500 

comparable  to  the  value  obtained  by  ref.  [28].  Meanwhile,  the 
methanol  feed  value  before  optimization  is  at  a  rate  of  1 .86 1  h-1 
while  after  optimization  is  at  the  rate  of  1 .42 1  h-1 .  This  shows  a 
good  improvement  of  the  system  after  the  optimization  process. 

Table  12  gives  the  results  of  the  system  efficiency  obtained 
in  this  study  as  compared  to  other  studies.  The  system  efficiency 
value  in  this  study  is  calculated  by  taking  into  account  the  fuel 
processing  system  efficiency  factor,  rjfp  as  0.95,  the  separation 
system  efficiency  factor  77 sp  as  0.95,  the  fuel  cell  stack  efficiency 
factor,  rife  as  0.65  and  the  parasitic  load  as  8%  of  the  output 
power.  The  system  efficiency  value  in  this  study  is  found  to  be 
within  the  range  obtained  by  other  studies.  Apart  from  that,  the 
system  efficiency  was  observed  to  have  increased  from  45.5% 
before  optimization  to  47.5%  after  optimization. 

The  next  paragraph  explains  on  model  verification  for  eco¬ 
nomic  parameters.  It  is  found  that  the  value  of  the  specific  stack 
cost  incurred  by  this  study  is  in  the  range  of  US$  500-700. 
This  value  is  comparable  to  the  ones  obtained  by  other  studies 
as  shown  in  Table  13.  However,  the  model  given  by  ref.  [45] 
expects  the  cost  to  be  as  low  as  US$  200-300  kW-1.  This  is 
because  the  cost  shown  is  referred  to  a  high  output  capacity 
of  up  to  500,000  vehicle  units.  Therefore,  it  can  be  deduced 
that  cost-saving  effort  is  very  much  dependent  on  the  output 
capacity.  This  has  been  put  forth  in  a  study  by  ref.  [46],  which 
predict  that  the  fuel  cell  cost  in  the  year  2020  will  be  at  a  rate  of 
US$  38  kW-1  with  an  annual  output  capacity  of  5  million  fuel 

Table  15 

The  comparison  of  specific  investment  cost  with  other  studies 


References  System  Cost  per  kW  (US$) 


[47] 

PEMFC 

2600-3000 

[44] 

PEMFC 

2500-3000 

[50] 

PAFC 

3000 

[51] 

PEMFC 

2500 

[29] 

SOFC 

2500-3000 

[52] 

PEMFC 

2000-3000 

This  study 

PEMFC 

2500-3000 

Table  16 

Infrastructural  investment  cost 


References 

Vehicles  capacity 

Infrastructural 

Vehicle  cost  (US$ 

production 

(million) 

investment  cost 
(billion  US$) 

per  kW) 

[47] 

25  (methanol-H2) 

45-95 

2000-3800 

This  study 

5  (methanol-H2) 

10-20 

2000-4000 

cell  vehicles.  The  estimated  specific  manufacturing  and  invest¬ 
ment  costs  obtained  by  other  studies  are  in  the  range  given  in 
Tables  14  and  15.  The  infrastructural  investment  cost  is  as  shown 
in  Table  16.  The  infrastructural  investment  cost  in  this  study  is 
in  the  range  of  US$  10-30  billion  with  an  output  capacity  of  5 
million  and  a  cost  of  US$  2000-4000  kW- 1 .  This  is  comparable 
to  the  study  by  ref.  [47]. 

7.  Conclusions 

The  main  objective  of  this  study,  which  was  to  solve  the  syn¬ 
thesis  and  optimize  the  design  for  a  PEM  fuel  cell  system  via  a 
separator-reactor  using  the  SQP  approach  for  non-linear  equa¬ 
tions,  was  successfully  obtained.  Although  the  non-linear  SQP 
approach  is  a  well-known  method  in  chemical  system  designs, 
this  study  is  the  first  attempt  for  a  fuel  cell  system.  The  complex 
equations  involved  have  been  simplified  into  matrix  forms  and 
solved  using  MATHLAB  software. 
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